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The empirical absence to date of particles obeying parastatistics in high energy collider exper¬ 
iments might be due to their large masses, weak scale couplings, and lack of gauge couplings. 
Paraparticles of order p = 2 must be pair produced, so the lightest such particles are absolutely 
stable and so are excellent candidates to be associated with dark matter and/or dark energy. If there 
is a portal to such particles, from a new scalar Ai boson they might be cascade emitted as a pair of 
para-Majorana neutrinos as in Ai —> A 2 i' a vp or as a pair of neutral spin-zero paraparticles such as in 
Ai —> A 2 AB, where B is the anti-paraparticle to A. In this paper, for an assumed supersymmetric- 
like “statistics portal” Lagrangian, the associated connected tree diagrams and their parastatistical 
factors are obtained for the case of order p = 2 parastatistics. These factors are compared with 
the corresponding statistical factors for the analogous emission of a non-degenerate or a 2-fold de¬ 
generate pair which obey normal statistics. This shows that diagrams, and diagrammatic thinking, 
can be used in perturbatively analyzing paraparticle processes. The parastatistical factor associated 
with each diagram does require explicit calculation. 


I. INTRODUCTION 

In the standard model all particles are either fermions 
or bosons which correspond to order p = 1 parastatis¬ 
tics. Identical fermions (bosons) occur only in the 1- 
dimensional totally antisymmetric column (totally sym¬ 
metric row) representations of the permutation group. 
Parastatistics is a natural and simple generalization 
which includes the additional higher dimensional repre¬ 
sentations of the permutation group. Fields and quanta 
obeying parastatistics are allowed in local relativistic 
quantum field theory [1-8]. Occasionally in this paper 
there are brief summaries, such as in the appendices, so 
as not to assume that the reader has a quantum field 
theory background in parastatistics. 

In this paper, we concentrate on order p = 2 parastatis¬ 
tics, which is the simplest such generalization of normal 
Fermi and Bose statistics. A simple consequence of or¬ 
der p parastatistics is that up to p identical parafermions 
(parabosons) can occupy a totally symmetric (antisym¬ 
metric) state, unlike for normal statistics. More gener¬ 
ally, identical parafermions (parabosons) of order p occur 
in Young diagrams with at most p columns (rows). 

Due to p = 2 parastatistics, an even number of para¬ 
particles must occur in the “total external state” for a 
physical process, so paraparticles must be pair produced 
and the lightest paraparticles are absolutely stable. The 
“total external state” consists of the particles in the ini¬ 
tial state plus the final state. Because of this absolute 
stability, paraparticles of order p = 2 are excellent can¬ 
didates to be associated with dark matter and/or dark 
energy (accelerated expansion), given what is currently 
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known from astrophysics and accelerator experiments. 

If there is a “statistics portal” from normal bosons 
and fermions to p = 2 paraparticles at a high energy 
collider, then these particles might be emitted in a cas¬ 
cade process from a new scalar Ai boson as a pair of 
para-Majorana neutrinos as in A± —> A 2 i' OL vp or as a 
pair of spin-zero paraparticles such as in A\ —> A 2 AB, 
where B is the anti-paraparticle to A. The paraparti- 
cles/parafields are denoted by a “breve” accent. All the 
new particles considered in this paper are assumed to 
be electromagnetically neutral with ~ 100 GeV to ~ 2 
TeV scale masses. The diagrammatic parastatistical fac¬ 
tors are calculated for these two pair emission cascades 
because of their massive and unstable final A 2 normal 
spin-zero boson, versus the empirical difficulties for inves¬ 
tigating a cascade to an almost massless final Majorana 
neutrino v 2 in A\ —> 1 / 2 AP. Depending on the unknown 
masses and coupling constants, these cascade processes 
might occur in the on-going experiments at the LHC with 
Vs ~ 13 - 14 TeV. 

As in the supersymmetric Wess-Zumino model [9], we 
assume that the portal Lagrangian densities for the cas¬ 
cade processes involve both a Majorana spin-1/2 field £ 
and a neutral complex spin-zero field A which respec¬ 
tively obey Fermi and Bose statistics, and parafermi and 
parabose counterparts £ and A which obey order p = 2 
parastatistics. We consider this complex A field in the 
particle-antiparticle basis with corresponding quanta A 
and B. Similarly, A and B are the quanta for the com¬ 
plex spin-zero parabose A field. We will assume that 
there are two new A\ 2 (with antiparticle B 12 ) bosons 
with nr 1 > m 2 >> 0, that all mass values are at the 
~ 100 GeV to ~ 2 TeV scale, and that each of the cascade 
processes is kinematically allowed. We also assume that 
if not for their weak-scale portal associated couplings, the 
paraparticles would only interact gravitationally. Obvi- 
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ously, the 7 cascade processes considered in this paper 
are kinematically analogous to —> prv„v T . However, 
the Hi^’s and Hi^’s are spin-zero, so there do not ex¬ 
ist useful polarization observables due to the cascading 
particle’s spin, but concurrently there are fewer unknown 
possible covariant couplings. 

Using these Lagrangian densities, we perturbatively 
calculate the S-matrix elements for A\ —> A 2 v a i>p, - ■ ■ 
and A 1 —> A 2 AB, • • •. We find that the tree diagrams 
for the associated connected amplitudes for cascade emis¬ 
sion of a pair of paraparticles correspond to the same dia¬ 
grams as in the case of the emitted pair obeying ordinary 
statistics, see Fig. 1 and others below. While the dia¬ 
grams are the usual covariant perturbative ones, with the 
initial state on the left and the final state on the right, in 
labeling the virtual lines by A or B , the displayed time¬ 
ordering has been assumed. The arrows on the particle 
A (antiparticle B) scalar boson lines are correspondingly 
forward (backward) in time. Since A is a complex field, 
upon a time reversal of a time-ordered virtual line, ex¬ 
change A and B label. Unlike the spin-zero Higgs bo¬ 
son which is its own antiparticle, the neutral A field has 
distinguishable particle-antiparticle quanta. This same 
time-ordering property holds for time-ordered A and B 
virtual lines associated with the A field. In the figures, 
vertices and lines associated with the paraparticles are 
drawn heavy or “dark.” There are also “dark dots” on 
the external paraparticle legs which enables omission in 
the figures of an awkward “breve” accent on the Weyl 
spinors. In the case of p = 2 parastatistics, the paras- 
tatistical factors c p for the diagrams displayed are evalu¬ 
ated. 

These c p factors in the p = 2 para case are then com¬ 
pared with the analogous statistical factors Cd calculated 
for the amplitudes in the case of the emitted neutral pair 
obeying ordinary statistics and in the case when there is a 
hidden 2-fold degeneracy, for instance A\ —> A 2 iy a ,a^a, 0 , 
where there are two kinds of emitted pairs v a ,aVa,fi with 
a = 1,2 the degeneracy index. In the 2-fold degenerate 
case, as for a final particle polarization summation, this 
index a is summed over to obtain the partial decay width. 
The assumed portal Lagrangian densities considered for 
these two comparison cases are analogous to those for the 
para case. 

In agreement with what might have been anticipated 
by some readers, our explicit calculations show that for 
each diagram the statistical factor c p for order p = 2 
parastatistics, and hence the associated partial decay 
width, is the same as the Cd statistical factor for such 
a 2-fold degeneracy. 

Section II contains the supersynnnetric-like La¬ 
grangian densities assumed for these cascade processes. 
It continues with the evaluations of the statistical factors 
c p in the para case and of the analogous factors Cd in the 
cases of emission of a non-degenerate or a 2-fold degener¬ 
ate pair obeying normal statistics. Section III discusses 
the predictions for partial decay widths for these three 
cases. Section IV has some concluding remarks. 



FIG. 1: Reading left-to-right, the first 2 of 4 diagrams for a 
cascade from a new scalar A\ boson by emission of a pair of 
para-Majorana neutrinos Ai —> Aii'cVp. The virtual scalar B 
is the antiparticle to A. Scalar bosons are denoted with thin 
dashed lines. In the diagrams in this paper, “dark dots” de¬ 
note the portal Lagrangian vertices and the external p = 2 
paraparticles which have weak-scale portal associated cou¬ 
plings. The “dark” solid lines for the para-Majorana neu¬ 
trinos are arrowed per the forward left-handed ad and back¬ 
ward right-handed y Weyl spinor, final state wave functions 
as in DHM [10]. The Greek subscripts label the momenta and 
helicities of the para-Majorana neutrinos in (17) and (18). 


The relatively simple tri-linear relations for the cre¬ 
ation and annihilation operators for an “order p = 2 
family” of parafields are listed in Appendix A. 


II. CASCADE PROCESSES WITH EMISSION 
OF A PAIR OF PARAPARTICLES 

A. Lagrangian densities 

For each of the interaction Lagrangian densities there 
is an explicit normalization of its coupling constant: For 
fields obeying normal statistics, a factor of (1/n!) occurs 
when that field occurs to the nth power. For a para La¬ 
grangian density, two parafields occur in their appropri¬ 
ate commutator/anticommutator ordering, see after (9), 
and also with an additional factor of (1/2). 

While these are the usual normalizations associated 
with the identity of the fields in normal statistics and 
in parastatistics, these definitions are arbitrary. How¬ 
ever, these definitions of coupling constants are fixed and 
are used to calculate the statistical factors ( c p and Cd ) 
for each diagram/amplitude. Any overall minus sign, or 
phase, is absorbed into the amplitude so c p ,Cd > 0. From 
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the values obtained for these factors, the consequences of 
alternate normalizations can be easily considered. The 
overall sign of each of the interaction Lagrangian densi¬ 
ties has been arbitrarily chosen as minus. 

Among the usual p = 1 fields, we consider inter¬ 
actions as in the supersymmetric Wess-Zumino model 
[9], but with unrelated weak scale coupling constants, 
so only slightly more general than in the supersymmet¬ 
ric limit. We use the excellent supersymmetric for¬ 
malism/notation of Dreiner-Haber-Martin (DHM) [10] 
with additional “breve” accents to denote the paraparti- 
cles/parafields. The fields have their usual covariant 
momentum-expansions and normalizations in terms of 
their associated creation and annihilation operators [11]. 
The interaction densities involving only p = 1 fields are 


Cy = - f -{A^ + A^U) 

(1) 

Cc = -\{A{A^) 2 + A 2 A^) 

(2) 

C q = ~{A*) 2 A 2 

( 3 ) 

For the cascade processes, we consider the following 
“statistics portal” couplings between these p=l fields and 
the p=2 fields, with anticommutator curly braces and 
commutator square brackets: 

Cy = - {([UlA + A%b 

( 4 ) 

C 2 c = -\{{A,A}A' + A{A\A'}) 

( 5 ) 

£ 3 c = -|(A + A t )M,A t } 

(6) 

C 2q = ~({A,A}(A') 2 +A 2 {A\A'}) (7) 

Ch = ~{AA'}{A,$} 

( 8 ) 

c a = ^{Ia} + {a\M 

( 9 ) 

In these Lagrangian densities, the 
rules of paraquantization dictate the 

standard 

connnuta- 


tor/anticommutator ordering of the parafields. By 
“paralocality” [3, 12] for fields obeying order p = 2 
parastatistics, two parafermi fields occur in a commuta¬ 
tor ordering, whereas two parabose fields, or a parabose 
and a parafermi field, occur in an anticommutator 


ordering. Paralocality is a generalization of locality for 
parafields, see Appendix B. 

For comparison, we also consider the case of cascade 
decays by pair emission fields A a (neutral complex spin- 
zero) and £ a (Majorana spin-1/2) obeying respectively 
Bose and Fermi statistics, for instance A\ —> A 2 iy a ,a^a, 0 - 
For the degenerate case, the Lagrangian densities are 
analogous to the above portal ones: 

£ d y = -^(AUa + A'U a ) ( 10 ) 

C d 2c = - t j(AU a A a +AlAiA) (11) 

C d 3c = -T d (AA a A{ + A a A{A^ (12) 

4 q = -^(A a A a (A^ + A 2 AiA\) (13) 

4 q = ~G d (A a Al)(AA') (14) 

£ d A = -jd(ttaA a + AU a 0 (15) 

For the case of 2-fold degeneracy, the degeneracy index 
(a = 1,2) is summed over in these densities. 

The interaction Lagrangian densities which do not oc¬ 
cur in the cascade processes calculated in this paper are 
(1) and (3), and in the paraparticle portal case (9) and its 
analog (15) in the degenerate case. However, the empir¬ 
ically difficult to observe cascade to an almost massless 
final Majorana neutrino i/ 2 in Ai —> v^Av does involve 
both (1) and the portal coupling (9), and its degenerate 
counterpart (15). 

B. Parastatistical factors for 7 cascade processes 

The above interaction Lagrangian densities have a 
particle-antiparticle transformation symmetry such that 
the results obtained for each cascade also hold for the cas¬ 
cade obtained by transforming all Ai •<-»• Bi and A o B. 
For instance, the parastatistical factors are the same for 
A\ —> A 2 AB and B\ —> B 2 BA. For the normal statis¬ 
tics cascades involving A a and £ a , there is the analogous 
transformation of all Ai o Bi and A a o B a . Conse¬ 
quently, the statistical factors c p and Cd obtained below 
for the diagrams in the scalar A\ decay process are the 
same as for the associated antiparticle B\ decay process 
because of this particle-antiparticle transformation sym¬ 
metry. 

1. Emission of a pair of para-Majorana neutrinos: 

A\ A^VaVp and A\ —> BiVaVp 

In this paper the evaluations of the S-matrix elements 
only involve processes with a pair of final paraparticles. 
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We calculate the associated amplitudes in the “occupa¬ 
tion number basis” for a specific ordering of the two parti¬ 
cles in the pair and then by addition or subtraction, con¬ 
struct the corresponding amplitudes in the “permutation 
group basis” [2] to obtain the physical amplitudes for the 
pair of paraparticles. In these evaluations, calculating in 
the occupation number basis halves the number of terms, 
versus using the permutation group basis, and a simple 
relabeling in the final expression gives the amplitude for 
the opposite ordering of the final two paraparticles. This 
distinction between fundamental bases in parastatistics 
is explicitly and simply explained below in the context 
of the calculation of the amplitudes for A\ —> A 2 i> a vp 
associated with the two Fig. 1 diagrams. This leads to 
the discussion in the text of the two physical permutation 
group basis final states of (19) below. 

In canonical quantum field theory, for particles obey¬ 
ing normal statistics there is a successful normal order¬ 
ing procedure for correctly ordered Lagrangian densities 
which is used in the perturbative evaluation of S-matrix 
elements [13]. This procedure discards various diagrams 
and yields results for the standard model which are cur¬ 
rently in highly precise agreement with experimental 
data. However, this procedure has not been generalized 
for paraparticles. Nevertheless, as shown in this paper, 
knowing from purely p = 1 quanta the canonical assem¬ 
bly of contributions from the perturbative evaluation into 
physical amplitudes, we find that it is straight-forward to 
proceed analogously by hand for p = 2 paraparticles us¬ 
ing the above paraquantized Lagrangian densities: We 
require each field in Li n t to contract with a field in a dif¬ 
ferent Li n t or with a particle in the initial or final states. 
This omits disconnected diagrams and ones with a single 
Li n t term self-contraction. In this context, it is important 
to note that there are highly non-trivial signs associated 
with this diagrammatic application of the tri-linear quan¬ 
tization relations and of the paralocal Lagrangian densi¬ 
ties involving order p = 2 fields. Clearly, the two crucial 
tests of this systematic diagrammatic evaluation of para- 
particle S-matrix elements will be whether it generalizes 
in perturbative quantum field theory and whether the 
resultant amplitudes do indeed agree with experiment. 

(i) We first consider a cascade from a new A\ boson 
by emission of a pair of para-Majorana neutrinos A\ —> 
A 2 v a vp: 

From and Lc there is the following time-ordered 
product 

Sfi = (*) 2 J d 4 x i J d 4 x 2 0(t\ - t 2 ) 

Aiv a vp\{L^(xi)Lc{x 2 ) + L c {xi)Lftx 2 )}\A k > 

(16) 

The final state has the z/ a z/g paraparticle operators or¬ 
dered as | Aiv a vp >a= > in the occupation 

number basis. The role of the subscript on the ket-state 
(bra-state) is to denote the place-position order [14]. We 
are labeling one place-position order “A” {a) ft) and its 


orthogonal counterpart “B” {ftft), where the ordering 
of the operators is reversed | A[V a vp >b= \ftftft\Q >. 

In p = 2 parastatistics there are the tri-linear relations 
instead of the usual bi-linear ones, so these A and B 
orderings in the occupation number basis must be distin¬ 
guished. 

We label the Ai creation operator by U. Notice the 
essential and easy to forget two paraparticles’ factor of 
(-^j) 2 in the state norm due to the vacuum condition 

a.kOt \|0 >= 2(5fcz|0 >. These extra paraparticle normal¬ 
ization factors occur because our calculations depend on 
the “arbitrary p normalization,” see Appendix A. 

By writing the fields of the A-ordered final state in 
(16) in terms of their positive- and negative- frequency 
parts, and then using the p = 2 tri-linear relations for 
the paraquanta, we obtain amplitudes corresponding to 
the two connected tree diagrams shown in Fig. 1. See 
Appendix C for p = 2 normalization details. 

From (16), the (si) amplitude [11] for the A-ordered 
final state is 

- iM^ 1} = {1 }{it)(if )~2 - 1 2 . y A (Pa, K)va{P 0 , \p) 

q A — + it 

(11 

with A being the 2-valued summed index for the com¬ 
muting two-component “right-handed Weyl” spinor, fi¬ 
nal state wave function tja of DHM [10], and the (s2) 
amplitude is 

- iM { A 2) = {l}{it){if)— - 1 . xftp a ,X a )x^ A {pp,Xft 

Qb — m + it A 

0 

with A-dot being the 2-valued summed index for the com¬ 
muting, conjugate “left-handed Weyl” spinor, final state 
wave function x^ A . Because of our usage of Greek let¬ 
ters for parafermions, such as in the tri-linear relations 
in Appendix A, we use undotted and dotted capital Ro¬ 
man letters for these 2-valued Weyl indices in place of 
the lowercase Greek letters in DHM. 

As in DHM, the arrows on the two-component spinor 
lines correspond to fields with undotted (dotted) indices 
flowing into (out of) any vertex. The direction of an 
arrow, versus a vertex, for either a spin-zero or a spin-1/2 
line is unchanged upon any time reordering of a displayed 
diagram. From the Lagrangian densities (1), (4), (9) and 
(10), (15), the arrows on the lines for the scalar fields A 
and A are also into (out of) any vertex per the common 
into (out of) direction of the two spinor lines. 

To maintain simplicity of the expressions for the ma¬ 
trix elements, we omit the associated mixing matrices 
between the mass eigenstates and the interaction eigen¬ 
states for the external A\ 2 bosons. In this paper, the 
amplitude/diagrammatic normalization is for a single A, 

£, or A in the virtual propagators. Also, while the stan¬ 
dard model lacks sufficient CP-violation for the observed 
baryon and lepton asymmetries of the universe, we omit 
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explicit CP-violation formalism and possible mixing of 
the A 12 with the B 12 bosons. 

For comparison, in the case with pair emission fields 
A a and £ a obeying the usual Bose and Fermi statis¬ 
tics, the same amplitudes for (si) and (s 2 ) are obtained 
for the process A\ —> A 2 ^ a ,aTa „8 with |^4/^ a ,ct^a,/3 >= 
Aa4/3j|0 > except in place of the parastatistical factor 
{1 }(</), there is instead a factor of {1 }(t/d), where the 
respective statistical factors c p and Cd are given in the 
curly braces. In writing these statistical factors times 
coupling constants, we omit each (i) associated with the 
iCi n t vertex. This is the comparison amplitude for all 
fields obeying ordinary statistics for the Lagrangian den¬ 
sities given in (10-15). In the 2-fold degenerate case 
where there are two kinds of emitted pairs v a ,otVa,@, cal¬ 
culation of the partial decay width requires a factor of 2 
due to summing over the two final degenerate channels. 

For the orthogonal B-ordered final state, the same 
amplitude for (si), and similarly for (s 2 ), is obtained 
but with an opposite overall sign in comparison to the 
A-ordered final state, so that the permutation group 
basis amplitudes and M.^ 2 ^ for the symmet- 

ric/antisymmetric final states 


1 

7t v 


\AlU a U/3 > sym,asym — /— (| A; V 0 Vfj > A A | AiU a Up > ) 


(19) 


are respectively zero and \/2 times those for the A- 
ordering. Hence, from the values of the statistical fac¬ 
tors c p and Cd , if these were the only two diagrams, upon 
summing over the two permutation basis final states for 



FIG. 3: The 3 diagrams for Ai cascading to £>2 by emission 
of a pair of para-Majorana neutrinos Ai —> B 2 V a i , i 3 . There is 
a para-Majorana mass insertion in the ( t ) and ( u ) diagrams. 


the decay process A\ —> A 2 V a vp the partial decay width 
would be twice that for the corresponding normal statis¬ 
tics process A\ —> A 2 V a ,aV a ,i 3 with a non-degenerate pair. 
However, the p = 2 partial decay width would be the 
same as that for the case of emission of two kinds of 
pairs t'a.aZ'a,/? due to summing over these two degenerate 
channels. 

For the Ai —> A 2 V a vp cascade, there is also a contribu¬ 
tion from (£j )) 2 which corresponds to the two diagrams 
in Fig. 2. Again, for each diagram, the B-ordering gives 
the same amplitude, but with opposite overall sign versus 
the A-ordering. Also, again for the A-ordering, the ex¬ 
pressions associated with the diagrams are proportional 
in the case of paraparticles and the p = 1 (normal statis¬ 
tics) case of non-degenerate Majorana neutrinos. The 
contribution of the ( u ) diagram is minus that of the ( t) 
diagram with a and /3 exchanged. In the para case, the 
( t) diagram has a factor of ( 1 }(/) 2 , and in the p = 1 
case there is a factor of {1 }{fd) 2 instead. In the evalua- 
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A 2 / A 2 



FIG. 4: First 3 of 6 diagrams for the cascade Ai -r A 2 AB by FIG. 5: The remaining 3 diagrams for Ai -r A 2 AB. 

emission of a particle-antiparticle pair of scalar paraparticles 
AB. B is the anti-paraparticle to A. These scalar paraparti¬ 
cles are denoted by “dark” dashed lines with forward (back- a para-Majorana mass insertion contribution. The am- 

ward) in time arrows for the particle A (antiparticle B). plitude for the (u) diagram is again minus that of the (t) 

diagram with a and /? exchanged. For the (t) diagram, 

in the para case there is a factor of {1}(/) 2 and corre- 
tion for the para case, there is a factor of 2 which arises spondingly in the p = 1 fermion case a factor of {l}(/d) 2 . 
from transforming the position space propagator vacuum 
expectation value to momentum space, see Appendix C. 

(ii) As shown in Fig. 3, there is a similar cascade from 2. Emission of a pair of scalar paraparticles: 

Ai to the antiparticle B 2 by the emission of a pair of Ai -»• A 2 AB, A i —>■ B 2 A 3 A 4 , • • • 

para-Majorana neutrinos, Ai —> 

For each diagram, for the A-ordeiing the para ampli- j n ^he remaining 5 cascade processes, A\ —> A 2 AB, 
tude is proportional to that obtained in the case of ordi- Ai A3A4, ■ ■ ■ , a pair of scalar paraparticles are 

nary fermion Majorana neutrinos. Also for each diagram, emitted. For each process, the obtained A-ordered am- 

the B-ordered expression is of opposite sign to that of the plitudes can again be considered in terms of its covari- 

A-ordering, so the permutation group basis amplitude is an ^ diagrams which are displayed in the figures. These 

again the asymmetric one. A-amplitudes in the para case are again proportional to 

From £y and Cc , for the A-ordering there is a single those in the non-degenerate case in which there is a scalar 

(s) diagram with a parastatistical factor of (l}(f/). For pair emitted. In the following, for each diagram the re¬ 
tire analogous p = 1 cascade A\ —> B 2 is ata is at p, there is a spective statistical factors c p and Cd are listed, 

factor of {1 }(tfd)- The contribution from (£y) 2 involves For each diagram the same amplitudes are obtained 
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,b 2 



FIG. 6 : First 3 of 6 diagrams for the cascade Ai —> .B 2 A 3 A 4 
by emission of a pair of scalar paraparticles A 3 A 4 . The sub¬ 
scripts on and A 4 are for momentum labeling which dis¬ 
tinguishes the two identical parabosons. 


for the A-ordered and B-ordered final states. Therefore, 
in the permutation group basis the associated symmetric 
final state has an amplitude of \[2 times that for the A- 
ordering, and the amplitude for the antisymmetric final 
state vanishes. For the first cascade A\ —> A 2 AB with 
emission of a particle-antiparticle pair of paraparticles, 
the symmetric/antisymmetric final states are 


\A 2 jAB 

^ sym,asym — -^(,\A 2 ,iAB >,4 ±\A 2 .jAB > b ) 


( 20 ) 


with A-ordering and B-ordering of the kets \ A 2 ^AB >, 4 = 
|ZUt£t|o > and \A 2 jAB > b = \l^B^A^\Q >. 

(i) Fig. 4 shows the first 3 diagrams for the cascade 

A\ —)■ A 2 AB. 

Fig. 5 shows the remaining 3 diagrams: 

From C 3 g, there is the ( q ) diagram with a factor of 
{!}(— G) versus {1}(— Gd)- The minuses occur here be- 


B- 


Ai 4 x 

---■>- < A ^A, 

S' 

(S) ^2c 


B- 


A X ' 

---■> - ' 

^ * -A 

J~>—K 

2c 

B- 


(t2) 


Ap 


— > — 


(u2) 


A 

B V- 

0 'v 

2c 


FIG. 7: The remaining 3 diagrams for Ai —> B 2 A 3 A 4 . 


cause we omit each (z) associated with the )X m t vertex. 
From Cc and £ 3 £, the (si) and (s2) diagrams each have 
a factor of {l}(fT) versus {1 }(tTd). From the (£ 3 g ) 2 con¬ 
tribution, the (t 1 ) and ( u ) diagrams each have a factor 
of {1}(T ) 2 versus {l}(Td) 2 . Interestingly, there is only a 
single diagram contribution from (£ 2 c) 2 - This (f2) dia¬ 
gram has a factor of (l}(f ) 2 versus {1 }{td) 2 - 

(ii) The analogous cascade from A\ to the antiparticle 
B 2 , Ai —> B 2 A 3 A a , has the 6 diagrams shown in Figs. 6 
and 7: 

From C 2 q, there is the ( q ) diagram with a factor of 
{ 1 }(— F) in the para case versus a factor of { 1 }(— Fd) 
in the boson case. From Cc and £ 2 g, the (s) diagram 
has a factor of (l}(ff) versus (l}(ttd)- As shown, the 
remaining four diagrams arise from £ 2 g and £ 3 g. They 
are (tl), («1), (f2), and (m2). Each has a factor of {l}(tT) 
versus {1 }(t d T d ). 

(iii) The cascade from Ai to A 2 by A\ —1 A 2 A 3 Ai has 
the diagrams shown in Fig. 8 : 

From Cc and £ 2 g, the (s) diagram has a factor of 
(l}(ft) versus {1 }(ttd)- From £ 2 g and C 3 g, the (t) and (u) 




£ 


£ 


>2 

y 

A X 

,a 2 

/ 

^ >l' 


B s»B 3 


(s) 

«B 4 



FIG. 8 : The 3 diagrams for the cascade A\ —¥ A2A3A4 by 
emission of a pair of scalar paraparticles A3A4. 

diagrams each has a factor of (l}(fT) versus {1 }{tdTd)- 

(iv) If instead there is emission of antiparticle pair 
-B 3 -B 4 via the cascade A\ -A A2B3B4, there are the dia¬ 
grams shown in Fig. 9: 

From Cc and £ 2 g, the (s) diagram has a factor of 
(l}(ft) versus {1 }{ttd)- From £ 24 and £ 3 ^, the ( t ) and ( u) 
diagrams each have a factor of {lKtT) versus {1 }{tdTd). 
The interaction vertices £ 2 c and £ 3 ^ in the (f) and ( u ) 
diagrams are exchanged in Fig. 9 for emission of B3B4 
versus those in Fig. 8 for emission of A3A4. 

(v) For the cascade A 1 -A B2A3B4 by emission of A3B4 
there are the diagrams in Fig. 10: 

From Cc and £ 3 ^, the (s) diagram has a factor of 
{l}(tT) versus {l}(tTd). From second order in £ 3 g, the 
(t) and (u) diagrams each have a factor of { 1 }(T ) 2 versus 
{ 1 }(/'./)•• 

As briefly explained in Appendix D, the same parasta- 
tistical factor c p (as above) is obtained for each diagram 
in the alternate p = 2 normalization of Green and Volkov 
[1] for the tri-linear commutation relations. To achieve 



FIG. 9: The 3 diagrams for the cascade A\ —> A2B3B4 by 
emission of an antiparticle pair of scalar paraparticles B3B4. 

these same c p values, there is a necessary rescaling of each 
of the portal coupling constants in (4-9) by g, —> 2 g t . 


III. COMPARISON OF PREDICTIONS FOR 3 
CASES 

To compare the partial widths in the 3 cases, in the 
Lagrangian densities we assume the corresponding cou¬ 
pling constants involved in the cascade are equal in the 
para case and in the two p = 1 cases of a non-degenerate 
or a 2-fold degenerate pair. 

For the assumed portal Lagrangian, the scalar pair 
emission mode A\ -A B 2 B 3 B 4 is forbidden through 
quadratic order in the Lagrangian densities. Conse¬ 
quently, when viewed inclusively, the 5 scalar pair cas¬ 
cades from the new A\ boson separate into 3 processes 
with A\ —> A 2 + X versus 2 processes with Ai —> B 2 + X 
because there is the A\ -a A 2 B 3 B 4 cascade. 
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B 2 



FIG. 10: The 3 diagrams for the cascade Ai —¥ B2A3B4 by 
emission of a particle-antiparticle pair of scalar paraparticles 
A3B4. 


The diagrams for these cascade processes with emis¬ 
sion of a pair of scalar paraparticles, or of a pair of scalar 
Bose particles, do have common values for all their re¬ 
spective c p and Cd statistical factors. This enables fac¬ 
torization of these common-valued c p and Cd into overall 
coefficients. When such a factorization occurs, the partial 
decay widths for the para case versus that for emission 
of a non-degenerate pair are related by 

dT p -2 = 2 |c p /c,i| 2 d,T p -i (21) 

where the 2 permutation group basis final states have 
been summed in the para case. From this expression, 
two times the partial width is predicted for all paraboson 
pair cascades versus the Bose case of emission of a non¬ 
degenerate pair obeying normal statistics. 

Similarly, the partial width for the para case can be 
compared with that for the case of emission of a 2-fold 


degenerate pair 

dT ' p —2 = \c p /c d \ dT deg.pair (22) 

The same partial width is predicted for all paraboson 
pair cascades versus emission of a 2-fold degenerate scalar 
boson pair. 

In the supersymmetric limit, the mass of the para- 
Majorana neutrino v would be the same as that for the 
scalar paraparticle A and its anti-paraparticle B , but 
in nature the paraparticle spin-1/2 and spin-zero masses 
might be different. This might enable kinematic separa¬ 
tion of a cascade process with emission of a pair of para- 
Majorana neutrinos from one with emission of a pair of 
scalar paraparticles. In the case of mass degeneracy of 
the v and A particles, generalizations of some of the tech¬ 
niques which exploit the neutrino spin in r~ —> ^“rai¬ 
ntight possibly be used to separate the v a vp cascades 
from the A 3 A 4 cascade. 

For the two para-Majorana neutrino cascades, A\ —► 
A 2 V a vp and A\ —1 B 2 V a vp, an overall factorization of 
c p and Cd is also possible, and the partial widths in the 
the para case are twice (the same as) the corresponding 
partial width for emission of a pair of Fermi Majorana 
neutrinos (2-fold degenerate Fermi Majorana neutrinos). 


IV. CONCLUDING REMARKS 

This paper is focused on showing that diagrams, and 
diagrammatic thinking, can be used in perturbatively 
analyzing p = 2 paraparticle processes for an assumed 
supersymmetric-like “statistics portal” Lagrangian. If 
there is a portal to such paraparticles at the LHC, they 
might be cascade emitted as a pair of para-Majorana neu¬ 
trinos as in A\ —> A 2 V a vp or as a pair of neutral spin- 
zero paraparticles such as in A\ —> A 2 AB. The asso¬ 
ciated connected tree diagrams and their parastatistical 
factors are obtained above for these 7 cascade processes, 
through quadratic order in the Lagrangian densities. For 
each diagram, these explicit calculations show that the 
statistical factor c p for order p = 2 parastatistics and the 
corresponding factor Cd for a non-degenerate or 2-fold de¬ 
generate pair which obeys normal statistics, satisfy the 
easy to remember c p = Cd = 1 relation. 

These results complement general quantum field the¬ 
ory results for arbitrary order p, including the general¬ 
ization of the spin-statistics theorem to “particles of half¬ 
integer spin obey parafermi statistics, while particles of 
integer spin obey parabose statistics” [15]. 

Certainly the systematic diagrammatic procedure used 
in this paper, which builds on the successful normal or¬ 
dering procedure for p = 1 fields, needs to be shown 
to generalize, especially to higher order non-tree dia¬ 
gram processes involving both p = 1,2 fields. However, 
from the herein calculations, it is noteworthy that the 
commutator ordering of two parafermi fields in the La¬ 
grangian terms (as dictated by paralocality for observ- 
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ables) is in agreement with the nontrivial respective ab¬ 
sence (presence) of coupling in the permutation group 
basis amplitudes for the final state \Aiv a vp > S ym,asym 
in A\ —> A 2 V a vp and A\ —> B 2 v a vp. This occurs dia¬ 
gram by diagram. Likewise, the anticommutator order¬ 
ing of two parabose fields in the Lagrangian terms is also 
in agreement with the absence (presence) of coupling, 
again diagram by diagram, in the permutation group ba¬ 
sis amplitudes for two final scalar parabosons in a totally 
antisymmetic (symmetric) final state in the 5 cascade 
processes, A\ A 2 AB , A± B 2 A 3 A 4 , 

While the permutation group basis is always physically 
required in constructing the associated physical ampli¬ 
tudes for all parabosons or all parafermions in the ex¬ 
ternal final (initial) states, the convenient usage of the 
occupation number basis in the calculations in this pa¬ 
per also generalizes to more than two final paraparticles: 

In the case of more than two parabosons, the cen¬ 
tral idea of only evaluating one occupation number ba¬ 
sis amplitude for each diagram works. For instance, for 
4 final parabosons of order p = 2, the totally symmet¬ 
ric final state which uses the totally symmetric bracket 
is 

—>= —7=(a^cV + + 

aVg Vg 

+ a)Ad)b^ + 6^d^c^a^)|0 > (23) 

The state has 6 independent orthogonal terms. However, 
for each diagram only one amplitude needs to be calcu¬ 
lated in the occupation number basis, for instance, the 
amplitude for the <A 6 ' A $ | 0 > term. The other ampli¬ 
tudes easily follow by permutations of the mode labels. 

In (23) and in the other state expressions in this sec¬ 
tion, the states are normalized, but with the factor 
for each paraparticle omitted. Note, the reordering rela¬ 
tions [16] of Appendix E must first be used to reduce the 
4! = 24 terms from the totally symmetric left-hand-side 
of (23), to 6 independent terms. There are 6 independent 
orthogonal terms because the sum of the dimensions of 
the three permutation group irreducible representations 
for 4 parabosons is 6. This totally symmetric permuta¬ 
tion group row representation is 1-dinrensional and has 
an eigenvalue of 6 for P SU m- This directly physical P SU m 
operator is the sum of the pair particle-exchange opera¬ 
tors, see (A6) for three parabosons in Appendix A. 

This single occupation number basis amplitude for the 
aW c t<ff|o > term then also suffices for construction of 
the permutation basis amplitude for each of the other two 
permutation irreducibles. For the L-shaped representa¬ 
tion with dimension 3 and an eigenvalue of 2 for P SU m , 
there is the eigenvector — ^a^d^c^)|0 >. Fi¬ 

nally, for the box-shaped representation with dimension 2 
and an eigenvalue of 0 for P sum , there is the similar eigen¬ 
vector i([a^&Md' i '+&Wdi'c'’] — [a^d^b^c 1 +d^a^ l c l b ^})\0 >. 

In summary, for 4 final parabosons for each diagram 
there is one occupation number basis amplitude which 
requires evaluation. By permuting the external particle 


mode labels, this single amplitude then gives by superpo¬ 
sition the three amplitudes which each correspond to the 
three distinct permutation group basis irreducible repre¬ 
sentations (the totally symmetric, the L-shaped, and the 
box-shaped). 

For 4 final parafermions the amplitude evaluation pro¬ 
cedure and the state decompositions are very similar with 
the permutation group basis irreducibles having the same 
dimension but opposite sign of P sum eigenvalues versus 4 
parabosons. Again, only one occupation basis amplitude 
needs to be evaluated for each diagram. Appendix F con¬ 
tains independent basis states and P S um eigenvalues for 
up to 4 parabosons (parafermions). 

Appendix A: Tri-Linear Relations for a “p = 2 
Family” of Parafields 

In the calculations of the cascade matrix elements, the 
following tri-linear relations [1] for a “p = 2 family” of 
parafields, A(y) and are used with parabose opera¬ 
tors denoted with Roman letters and parafermi operators 
denoted with Greek letters. In the supersynnnetric-like 
model in the present paper, there are of course an equal 
number of parabose and parafermi degrees of freedom. 
The parastatistics term “p = 2 family” means that all 
the fields in the family mutually obey these tri-linear re¬ 
lations [7]. The fields, £(x) and A(y), have their usual 
covariant momentum-expansions and normalizations in 
terms of these creation and annihilation operators, see 
(Cl) and (C2) below. In the arbitrary p order tri-linear 
relations, versus the following p = 2 tri-linear relations, 
there are twice as many terms on the left-hand side of 
each relation due to an additional overall commutator 
ordering [3]. 

The mode index /c,/,to includes the momentum com¬ 
ponents, and the helicity components for the para- 
Majorana field A and the A, B particle-antiparticle dis¬ 
tinction for the A complex field. For instance, in the tri- 
linear relations below for the para-Majorana operators, 
the generalized Kronecker delta is Si m = d\ l x m S^{pi — 
p m )- Here, for clarity, we omit/suppress a possible but 
awkward “breve” accent which might be put on top of 
each of the creation and annihilation operators. 

Several simple patterns are apparent: As for the usual 
p = 1 bi-linear relations, in each relation the left-hand- 
side has the second term with the three operators written 
in opposite cyclic-order to that of the first term. The sec¬ 
ond term has a plus (minus) sign when mostly parabosons 
(parafermions) occur in the tri-linear relation. On the 
right-hand-side, the existence of a Kronecker delta term, 
and its sign, corresponds to an aka] or aka] adjacent-pair 
factor from the left-hand-side. The tri-linear relations 
maintain the associated odd (even) “place positions” [14] 
of both the mode and also of the parafermi/parabose la¬ 
beling of the operators, whether reading left-to-right, or 
right-to-left. These simple properties also occur in the 
adjointed relations. The normalization of these p = 2 re- 
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lations corresponds to that of the tri-linear relations for 
arbitrary p parastatistics [1, 3]. The usual p = 1 cre¬ 
ation and annihilation operators for boson fields, such as 
the scalar complex field A, commute with these p = 2 
operators and those for fermion fields, such as the Majo- 
rana spin-1/2 field £, commute (anticonnnute) with the 
parabosons (parafermions). 

For all parabosons (Roman letters): 

UkdlUm d m dld k = 0; 
akaia) m - a' m aiak — 26 i m cik 

CLk^ &rn d m di dk — 2 SklCLm 2 5 m ld k (Al) 

For all parafermions (Greek letters): 

OikOLlOLm H - — 0? 

akOLia' m + a^aiat = 2Si m a k 

OikOL^OLm T dimdi^Oik — 28klOZm T 2 S m lCXk (A2) 

For two parabosons and one parafermion: 

akdipm - /3 m aia k = 0, 
dkfildm d m [3ldk = 0, 
akdlPln - PjndlClk = 0 , 
dkpialn - alnfilCLk = 0 , 
fikdld^m djridlfik — 2 Slmftk 
d k dj (3 m fimd\dk — 28klfim 

dkfi\ dm dmfij CLk — 0 (A3) 

For two parafermions and one paraboson: 

-(- b m aia k — 0, 

QkblCXm “t“ O^mblCXk — 0, 
akaib^ + bmaiak = 0, 
akbia'm + a\ % biak = 0 , 
bkOaaln + a'maibk = 28 i m bk 
d^kd^i bm T bm^i &k 28klb m 

atkb\ ctm + a m b\ a k = 0 (A4) 

In this arbitrary p order normalization, the important 
associated vacuum conditions for any mode indices fc, l 
are 

a fe |0 >= a /|0 >= 0 , < 0|0 >= 1 
a k a\ |0 >= a k a] |0 >= 24/10 > (A5) 

and afcCt/|0 >= afca||0 >= 0. The parabose and 
parafermi number mode operators are respectively N k = 
\{ a \, a k} ~ 1 and Af k = \ [a\,a k ] + 1. 

For p order parastatistics, in the vacuum conditions 
(A5), there is the substitution 2 —> p, and in the number 
operators the ending terms =Fl —> =F§, so the important 
zero point energies scale with the order p. Associated 
with (A5), for p = 2 there is an extra (-^) factor for 


each paraparticle in an external state. Thereby, the scat¬ 
tering matrix, and associated in-going and out-going par¬ 
ticle fluxes, have a common “particle density per unit vol¬ 
ume” normalizaton [11] for all external particles whether 
of order p = 2 or p = 1. 

For an initial state, final state, or observable expressed 
as a function of creation and annihilation operators, the 
directly physical “particle permutations” are products of 
the pair particle-exchange operators j = P 3 , t which 
exchange the i and j identical particles, so oq O dj or 
a\ O a}y As in [14], in the present paper these operators 
are denoted with an “overbar.” Instead, “place permuta¬ 
tions” are products of the pair place-exchange operators 
P rs = P s r which exchange the occupants of positions r 
and s in a creation and annihilation operator expression 
regardless of the identity of the occupants. 

In p = 2 parastatistics, unlike for p = 1 quanta, the 
external “permutation group basis” states are in gen¬ 
eral not eigenstates of the pair particle-exchange oper¬ 
ators P t ,j- Indeed, for two identical paraparticles, the 
external states are pair particle-exchange eigenstates as 
in (19) for parafermions and in (20) for parabosons. 
For three identical parabosons, the totally symmetric 1- 
dimensional external state is also even under each of the 
three particle-exchanges P t j- However, the three identi¬ 
cal paraboson state corresponding to the 2-dimensional 
L-shaped representation has basis vectors which are not 
eigenstates of the three P t j- For this mixed represen¬ 
tation and using the reordering relations of Appendix 
E, this is apparent because the two independent basis 
vectors can be chosen as -^(a^b^c^ — 6^c^a^)|0 > and 

-^(a'f&i'c'i'-bfri’cW — 2c'i'a' f &' i ')|0 >. When one of the three 

Pij acts on either one of these two basis vectors, it gives a 
linear combination of them. More generally, acting with 
any of the Pi.j on a state in an n particle irreducible 
representation of the permutation group preserves its ir¬ 
reducible representation. 

The sum of the three particle-exchanges Pj j which we 
denote P sum 

P sum = Pa, b + Pb,c + P c,a (A6) 

lias respective eigenvalues 3 and 0 for these two parabose 
representations, (row) and (L-shape), and so it can also 
be used to label them. Similarly, we find that states in n 
dimensional parabose representations are eigenstates of 
Psum (at least thru the four 6 paraboson irreducibles). 
For the ?i-dimensional totally symmetric representations, 
the eigenvalue n(n — l)/2 is equal to the number of 
pair particle-exchange operators. For states of identi¬ 
cal parafermions, the eigenvalues of P S um for the cor¬ 
responding irreducible representations are negative. A 
diagonal mirror reflection of rows and columns trans¬ 
forms a paraboson irreducible representation to a cor¬ 
responding parafermi irreducible representation. While 
the dimension of the permutation group representation 
provides one label for the irreducible representation for 
p = 2 external state, the 2-particle state shows that the 
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Psum eigenvalue is also required for a unique labeling. 
Also for the 6 paraboson state, the P S um eigenvalues of 
9 (L-shape) and 3 (box-shape) of distinguish these rep¬ 
resentations which are both 5-dimensional. 


Appendix B: Paralocality, Green Components, and 
Possible Additional Interaction Terms 

At the beginning of Section 2, in the construction of the 
supersymmetric-like portal Lagrangian densities, “par¬ 
alocality” is used. It is a generalization of locality for 
parafields [3]. The allowed forms of paraparticle cou¬ 
plings arise as a consequence of the tri-linear commuta¬ 
tion relations and the assumed locality condition. By 
locality, for the two obserables 0 (x) and O (y), their 
commutator must vanish 

[O{x),o\y)]=0, x~y (Bl) 

when points x and y are spacelike separated (denoted 
by the symbol ~). In the interaction picture, 0(x ) and 
O (y) are polynomial functions of the free parafield op¬ 
erators (f>i(x) which act on the vacuum of the physical 
Hilbert space. In arbitrary order p parastatistics, “par¬ 
alocality” holds when (Bl) is valid in the larger Hilbert 

space of the Green component fields <j>\ a \x) defined by 
the expansion 


by a Ci (the hat accent denotes Ci or cj ) is 

:[ci,c 2 ,---c„] t : = n! ^ c^ ai) 4“ 2) ■ ■ ■ c ^ n) 

0*1 i a n 

(B3) 

where the summation is over all different values of the 
Green indices Oi, a 2 , • • • , a n . As denoted by its redun¬ 
dant =F subscript, this “dotted bracket” or “second unit 
observable” is totally antisymmetric (symmetric) with 
respect to the labels 1,2, ...n in the case of all para- 
bose (parafermi) operators and it respectively vanishes 
for n > p, which is another meaning for the order p. In 
permutations for (B3), it is understood that the dagger, 
or no-dagger, on Cj, moves with the subscript i. 

In the context of deriving the most general selection 
rules for particles obeying parastatistics, these second 
unit observables were introduced in [3], (see earlier [17]). 
Such additional terms are treated in detail in Ref. [7]. 
However, these additional terms are forbidden if either 
there is the stronger locality condition 

[4>i(x),0 (y)] = 0, x ~ y (B4) 

or if there is a global symmetry such that the Green in¬ 
dices transform under 0(2) or U(2), instead of the smaller 
SO(2) or SU(2). These two locality conditions, (Bl) and 
(B4), are equivalent for ordinary fermions and bosons but 
are not for even-valued orders of p, see [7]. 


Mx) = ^2^ a \x). (B2) 

a —1 

where a is the Green index. For parabosons 
(parafermions) these Green component fields with the 

same Green index, <ftf\x ) and ^-“^(p), obey the usual 
Bose (Fermi) commutation relations, but anticommute 
(commute) with all <jtf\z) and (z ) for a ^ b. In an 
“order p family” of parafields, a parabose Green compo¬ 
nent and a parafermi Green component have the same 
commutation pattern as two paraboson Green compo¬ 
nents. Green components were introduced in [1], In [12, 
7], it is shown that locality implies paralocality. 

While in the perturbative calculations in this paper for 
order p = 2 we do not expand in Green components, they 
are very convenient tools for analysis and for checking. 
Historically, Green components have been exceptionally 
useful in developing and understanding fields and quanta 
obeying parastatistics, especially for arbitrary p order. 
Their underlying presence in parastatistics is a strong 
physics/mathematics motivation for the consideration, 
throughout the present paper, of the comparison with 
cascade emission of a 2-fold degenerate pair of particles 
which obey normal statistics. 

In above portal Lagrangian densities, we do not con¬ 
sider possible additional “second unit observables” which 
are allowed by paralocality. A generic example, in terms 
of paraparticle creation or annihilation operators denoted 


Appendix C: Evaluation of p = 2 Matrix Elements 


Some care is needed in the evaluation of p = 2 ma¬ 
trix elements because of the factor of 2 in the vacuum 
condition Ckc\ |0 >= 2<5fc/|0 >. In the “arbitrary p nor¬ 
malization” of [3] which is used in this paper, a factor of 
does not occur in a parafield’s momentum-expansion 

such as for the complex spin-zero parabose A field 


A(x) = 


d 3 q 


\J ( 27 r) 3 2 w, 


--(A(q)e~ 


B\q)e iqx ) 


(Cl) 


Consequently, A(a;)Al(y)|0 >= 2S^(x — p)|0 >, and 
2iAp(q) = 2i/(q 2 — mr + ie) in momentum space corre¬ 
sponds to < 0|T(A(a/)Al(:r))|0 >, with these three 2’s 
replaced by p’s for p order parastatistics. Similarly, there 
are corresponding factors of 2 occurring for the Majorana 
spin-1/2 parafermi £ held 


U(x)= l JL =( x a(p, X)a\(p)e lp '• 

1 y/(2lT) d 2L0 p 

+y A {p,X)a\{p)e lp ' x 


A=±l /2 ' 


(C2) 


with £a( x ) ~ ?!i( a: ) = (£a(®))^- For essential proper¬ 
ties of the commuting two-component Weyl spinor wave 
functions xa and yA see DHM [10]. 
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Appendix D: Same Parastatistical Factors for the 
p = 2 Normalization of Green and Volkov 

A simple, but very partial, working check of the above 
perturbative evaluations is to use the alternate p = 2 
normalization of Green and Volkov [1] for the tri-linear 
relations of Appendix A. If it is used, the same parasta¬ 
tistical factor c p is obtained for each diagram for these 
7 cascade processes. The differences are: (i) For both 
the p = 2 parabosons and parafermions, the quanta op¬ 
erators Cj —> \/ 2 di in the tri-linear relations and in the 
vacuum conditions (A5) of Appendix A. (ii) Each of the 
portal Lagrangian densities is quadratic in the parafields, 
so each of the coupling constants in (4-9) must be rescaled 
by gi —> 2g t . The momentum-expansions in Appendix C 
are still used without additional factors of y/ 2 , so there 
is no change in the normalization of the parafields versus 
the arbitrary p normalization used in this paper. Con¬ 
versely, if this Green-Volkov p = 2 normalization is used, 
then the arbitrary p normalization is obtained by the sub¬ 
stitution di but with the extra -^= factor, which 

would occur in each parafield momentum-expansion, in¬ 
stead moved out to be with the coupling constant in the 
Lagrangian density. 

Appendix E: Reordering Relations for p = 2 

With upper (lower) signs for parabose (parafermi) cre¬ 
ation operators, the reordering relations [16] are 

a^b^c 1 = ±c^V a* 

(aWXcW) = (c t d t )(a t 6 t ) = 

±(c t 6 t )(a t d t ) = ±(a t d t )(c t 6 t ) (El) 

The 3 operator relation is a cyclic one. In the 4 oper¬ 
ator relations, the unnecessary pairing parentheses are 
for displaying the pairing patterns. These patterns are 
that in reordering: (i) the even (odd) place positions are 
maintained, whether reading from the left or right, and 
(ii) each single exchange of a* fA c* or 6' -o- d} gives a ± 
sign. Both patterns are also in the 3 operator relation. 
The annihilation operators satisfy the same relations (re¬ 
move daggers). 

As discussed in the last several paragraphs of the text, 
these relations are particularly useful in the construction 
of the independent orthogonal terms needed in the exter¬ 
nal states in the permutation group basis. 

These reordering relations follow from the first lines of 
(Al) and (A2). For 3 mixed parabose and parafermi, all 
creation (annihilation) operators there are also reorder¬ 
ing relations corresponding to (El) which follow from the 
first two lines of (A3) and of (A4): These analogous 3 op¬ 
erator relations to (El) are also cyclic and hold with the 
± sign for mostly parabosons (parafermions). The ap¬ 
propriate signs in the 4 or more operator relations for 
any mixture of parabosons and parafermions then follow 
iteratively. 


Appendix F: Independent Basis States and P sum 
Eigenvalues for Up to 4 Parabosons (Parafermions) 

In this appendix “eigenvalue” for an n parabose or 
parafermi state means the eigenvalue of the P SU m oper¬ 
ator which is the sum of the pair particle-exchange P t3 
operators, see end of Appendix A. As in the concluding 
Sec.IV, in this appendix we are suppressing the extra 
normalization factor for each paraparticle. When these 
extra factors are included, each state is properly nor¬ 
malized. Each term in these states is independent and 
orthogonal because the reordering relations of Appendix 
E have already been used. 

For 2 parabosons, both of the two irreducible represen¬ 
tations are 1-dimensional. These correspond to the to¬ 
tally symmetric row (antisymmetric column) in a Young 
diagram, and have basis states 

-^[a t ,6 t ]±|0> (FI) 

with eigenvalues ±1. The bracket’s ± subscript denotes 
anticommutator (commutator). 

For 3 parabosons, there is the totally symmetric 1- 
dimensional row representation with 

-|=(a^?>t c t-|-5t c t a t c t a tftt)|o > (F2) 

v3 

and eigenvalue 3. For the 2-dimensional, L-shaped per¬ 
mutation group representation with an eigenvalue equal 
to 0, two basis vectors are 

^(a t 6 t c t -6 t c t « t )|0> 
v 2 

—^(a^c* + — 2cW&^)|0 > (F3) 

V6 

For 4 parabosons, there is the totally symmetric 1- 
dimensional row representation with 

—i=(a^ b^ c 1 dt + + 

a^d^b^c 1 + + &^ c t a t)|o > (F4) 

and eigenvalue 6. For the 3-dimensional, L-shaped rep¬ 
resentation with an eigenvalue equal to 2, three basis 
vectors are 

-^=( 0 )c*d) — &Wdt c t)|o > 

— ^—{a^d^b^c^ — dWc*&t)|0 > 
v 2 

—^(a^c^&t — fet^t c t a t)|o > (F5) 

v 2 

For the 2-dinrensional, box-shaped representation with 
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an eigenvalue equal to 0, two basis vectors are 

i([a^ ftftft + ft ft ft ft] — [ft ft ft ft + dWcV])|o > 
-^=([ftftftft + ft ft ft ft] + [cftcftftft + (ft ft ft ft] + 

—2[ftftftft+ftftftft])\0> (F6) 

For the n parafermion states, first recall that a diag¬ 
onal mirror reflection of rows and columns transforms a 
paraboson Young diagram irreducible representation to 
the corresponding parafermi representation. The P SU m 
eigenvalue is minus that for the corresponding paraboson 
representation. For the 2 parafermion irreducible repre¬ 


sentations, there are the (FI) basis states but with Greek 
letters per the notation in Appendix A. 

For 3 and 4 parafermions, the basis vectors are the 
same as above but with a change from Roman to Greek 
letters. In the above basis vectors the expected different 
signs for parafermions, versus parabosons, have already 
been absorbed by using the reordering relations of Ap¬ 
pendix E to reduce the basis vector expressions to the 
displayed independent orthogonal terms. For instance, 
if for 3 parafermions, one constructs the 1-dinrensional 
totally antisymmetric representation by first explicitly 
writing out the 3! terms, the re-ordering relations can 
be used to reduce it to (F2), but with Greek letters. 
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